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Abstract 

Large-scale astrophysical systems are non-extensive due to their long-range 
force of gravity. Here we show an approach toward the statistical mechan- 
ics of such self-gravitating systems (SGS). This is a generalization of the 
standard statistical mechanics based on the new definition of entropy; 
Tsallis statistical mechanics. Developing the composition of entropy and 
the generalized Euler relation, we investigate the galaxy distributions in 
count-in-cell method. This is applied to the data of CfA II South redshift 
survey. 



1 Introduction 

Astrophysical systems in the Universe are characterized by the gravitation. 
The structure formed through this long-range force is quite different from those 
formed through the other short-range forces. If the system does not strongly 
depend on the initial conditions of the Universe, we can apply statistical me- 
chanics for describing such self-gravitating systems (SGS). However, we cannot 
directly apply the standard Boltzmann statistical mechanics for SGS since the 
long-range nature of gravity strongly violates the extensive property of the sys- 
tem which is the premise of statistical mechanics. Actually, the total energy 



1 



increases much faster than the particle number N, the partition function Z 
often becomes complex]!]], reflecting the fact that there is no absolute stable 
state in SGS. 

In order to seek for workable statistical mechanics of SGS, we try an ap- 
proach based on the new definition of entropy whose extensivity is violated from 
the beginning; Tsallis statistical mechanics. 

We formulate the count-in-cell method for the large scale galaxy distribu- 
tions in this new statistical mechanics. First we calculate the expression of 
the composite entropy and the generalized Euler relation in this new statistical 
mechanics. These are applied to the data of CfA II South redshift survey. The 
parameter q becomes negative, which represents the instability of gravity. 

2 Tsallis Statistical Mechanics 

The ordinary Boltzmann statistical mechanics is characterized by the entropy 
S = —J2Pi^ n Pi- The distribution function pn,e = exp [— (E — /j,N) /T] /H 

i 

maximizes this entropy with constraints of the probability conservation, the 
energy conservation, and the particle number conservation. This statistical 
mechanics is originally aimed to describe the multi-fractal and chaos structures. 

It is characterized by the entropy of the form Q S q [p] = \J2pj — 1 J /(I — q), 

where q is a real parameter. Tsallis distribution function is obtained so that it 
maximizes this entropy with the same constraints. The solution has a power 
law tail: 

PN>E = ±{l-^-(E-E-»(N-N))}~\ (1) 
where T = CT, and C = J2p1- Partition function is defined as 

i 

Z^^ll-^iE-E-^N-N))}^ 9 . (2) 

N,E 1 1 J 

Note that pn,e reduces to the ordinary Boltzmann form for q — ► 1. For the 
consistent formulation, it is important to notice that the observable expectation 
value is calculated by the escort distribution ||: Pj = pj/C, (Q) = PiQi- 

i 

The averaged quantities E, N in Eqs.(|l]) and @) should be understood in this 
sense. 

3 Galaxy Distribution in Count-in-Cell Method 

There are many approaches to describe large scale structure of the Universe. 
In this paper, we will concentrate on distribution of galaxies. One of the well- 
known method to describe the distribution of galaxies is the two-point correla- 
tion function. However, when we solve equations of motions for the two-point 
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correlation function, we need three-point correlation function. Such higher or- 
der essentialness is called BBGKY chain. Since we have to cut the BBGKY 
chain, we have to apply some kind of approximations. 

On the other hand, count-in-cell method is often used to describe distri- 
bution of galaxies. In this method, we use an analytic formula for probability 
f(N) of finding N galaxies in a randomly positioned volume V. As we don't 
have to use any approximation in count-in-cell method, we can apply it even 
for clustered system. Saslaw and Hamilton [Q], and S. Inagaki introduced the 
virial parameter b= (gravitational correlation energy) /(kinetic energy of random 
motion), which measures the deviation from the dynamical-equilibrium. Then 
they found that in thermal-equilibrium, their theoretical investigation of f(N) 
is consistent with he N-body simulations or catalogues of observations. Strictly 
speaking, we should not apply thermal-equilibrium theory for expanding Uni- 
verse. However their consistency let us further study thermal-equilibrium sta- 
tistical description. In evolution of the Universe, dynamical-equilibrium 6=1 
must be realized before the Universe reach thermal-equilibrium. Therefore in 
this paper we consider the dynamical-equilibrium 6 = 1 case. We believe that 
there exist adequate fitting parameter other than the virial parameter 6. That 
is the reason why we consider non-extensive statistics which gives us parameter 

q- 

Supposing the galaxies distribute according to Tsallis statistics, we consider 
a system described by equilibrium thermodynamics. That is the grand canonical 
ensemble, characterized by the given temperature and the chemical potential. 

The probability to find no galaxy in the volume V is 

/(0)^E^ = J Po,o = %^, (3) 



where 



Po,o — 1 



''• • (4) 



The partition function can be decomposed as 



Y, (l-^J-(E-E-»(N-N))) =fe 9 ) c (5) 



q 

N,E 

and therefore E q = CVC 1 ^) = (1 + (1 - q) S) l/(1 ~ q) . Thus we obtain 

/ (0) = (1 + (1 - q) S)& (l + (E - /iff)) ~ q . (6) 

In order to reduce the above expression, we need to calculate the composite 
entropy and the Euler relation in Tsallis statistics. When we compose two 
systems A and B, the distribution function is given by p i j(A,B) = p i {A)pAB) 
and the composed entropy is S A+B = S A + S B + (1 — q) S A S B . Sequentially 
using this composition S N+1 = s + S N + (1 — q) sS N , (where s means entropy 
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of 1 system), we obtain the total entropy of N identical systems as 

(l + (l- g ) a )"-l 
*N ~ — q • (7) 

Generalized Euler relation is given by the following arguments. The vari- 
ables E,V,N are the natural arguments of the entropy: S N = S (E,V, N). 
Differenciating S aN = S (aE, aV, aN) with respect to a and setting a — ► 1, we 
obtain 



dS aN 



da 



OS (aE, aV, aN) 



i T T T ' 



+1 da 

where we used dS/dV = p/T, dS/dE = 1/T, dS/dN = -fi/T, which are 
guaranteed by the Legendre structure of the Tsallis thermodynamics ||. In the 
above, we can put any other values for a as well. In general, the a-dependence 
is inherited from S to T on the right hand side of Eq.@. Non-extensivity of S 
necessarily accompanies the non-intensivity of the Legendre-conjugate variable 
T. Actually, the fact that the variables E, V, N are extensive and p, fi are 
intensive in Eq.(^) renders the a-dependence of the temperature: 

T(a)=T 1 (l + (l-q)s) l - aN (9) 

where 7\ is the temperature of one particle system. Note that T = CT defined 
just after Eq.([l]) is a- independent! Probably this quantity T should be related 
with the physical temperature as defined from the velocity dispersion of the 
system. However at present, we do not have idea on the meaning of the quantity 
T0. This point will be further discussed in the last section. 

Then we obtain the Euler relation in non-extensive statistical mechanics 

(l + (l-g)S)ln[l + (l-g)5] = E + pV-^N 

1-q T 1 ' 

This temperature on the right hand side is T (a = 1) in Eq.(|9|), though the 
explicit form of the temperature does not appear in the final expression of 
f (N). Using this Euler relation, we can now eliminate E — fiN in / (0) and 
Eq.(||) reduces to 

q 

f (0) = (1 + (1 - q) (l - C (1 - q) + In [(1 + (1 - q) 5)]) ^ . (11) 

Since we consider the dynamical-equilibrium system, the system is fully virial- 
ized 6 = 1, and therefore the pressure p must be 0, we obtain: 

/ (0) = (1 + (1 - q) S)T^ (1 + In [(1 + (l-q) S)])^ . (12) 
Moreover using Eq.(R), we finally obtain 

-AT q 

/(0) = {1 + (1 - q) s}— [1 + N In {1 + (1 - q) s}}~ . (13) 

1 This situation is similar to the argument on the physical probability distributions at the 
end of the section two. Among P 4 and p t , related with each other by the factor C, we have 
chosen the escort distribution P i as the physical distribution. 
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Our parameters are s (the unit of non-extensive entropy per galaxy) and q (the 
Tsallis statistical parameter). 

Probability of finding N galaxies in the volume V is generated from / (0):^ 

(- n ) N d N 

where n is the galaxy number density ||. This is because the void probability 
/ (0) contains all the information of the whole correlation functions: 

f(0) = l-jp{X 1 } + j JpiX.Xzj-J J J 'P{X 1 X 2 X 3 }±..., (15) 

where P{X 1 X 2 } is the probability that there are galaxies in dVi at X l and 
dV 2 at X 2 . The above expression guarantees that the probability is properly 

oo 

normalized: f C^O = 1- 

N=0 



4 Comparison with observations 

We use the data of CfA II South redshift observations which includes 4392 galax- 
ies M. We have to reduce the data to the uniform sample. First we restrict the 
data to the galaxies whose absolute luminosity is brighter than the magnitude 
-19.1 and the distance 4,000 ~ 8, 000 [km/ sec]. The distance is measured by 
the cosmic redshift. We further exclude the edge of the observation region. We 
applied the K-correction for compensating the reddening. Finally the data is 
reduced to 870 galaxies. The number density is n = AAA x 10~ 9 [(km/secr 3 ]. 

We first fit the void probability / (0) by varying the parameters q and s. 
The best fit is realized by q = -5.66847 and s = 0.164142 (FigJT|). We fix these 
values and do not change them hereafter in this paper. With these parame- 
ters, general probability / (N) is given. In Figs.|]-||, we compared our calcula- 
tions and the CfA data. We have checked the normalization of probabilities 

oo 

J2 f (N) = 1 is realized. 

The negative value of q we obtained may not be so surprising. In refer- 
ence H, the value of q appears to decrease monotonically from 1 to — oo for 
one-dimensional logistic maps. Their model reveals unstable onset-to-chaos 
attractor. For negative q, the entropy functional loses its convexity and the dis- 
tribution becomes unstable. Thus the intrinsic instability of SGS is faithfully 
represented in this formalism and this fit. 



5 Conclusions and Discussions 

We have constructed the non-extensive statistical mechanics based on the non- 
extensive entropy. Especially calculating the entropy of composite systems and 
deriving the generalized Euler relation in thermodynamics, we could evaluate 
the void probability function / (0) and the probability to find N galaxies / (N). 

2 Note that this expression is slightly different from that by Saslow and Hamilton Q. 
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Figure 1: The void probability /(0). CfA II South Observations (dots with 
error bars) and our calculations (Solid lines) are plotted. From this fit, we 
obtained the parameters q and s, which should be kept unchanged hereafter. 



This result was applied to the CfAII South galaxy observations and we have 
obtained negative parameter q. This is thought to be another representation of 
the intrinsic instability of SGS. It will be also interesting to notice the fact that 
the multi-fractal scaling is observed in this CfAII South data within the scale- 
region from 500 [km/sec] to 3000[km/sec]||. We would like to clarify possible 
connection between the non-extensive distributions and the multi-fractal nature 
in the context of gravity. 

On the way we derive / (N), we encountered "scale (a) dependent temper- 
ature T". If we put the values we obtained q = —5.66847 and s = 0.164142 to 
Eq.(^), we can explicitly plot the scale dependence of the temperature. It turns 
out to reduce with increasing scale and abruptly drops to zero at about r rs 
600[km/sec] or, assuming the cosmic expansion speed as H = 72[km/sec/Mpc], 
at about 8.3[Mpc], which is almost the scale that the galaxy correlation func- 
tion becomes unity. On the other hand, it is apparent that the galaxies do have 
peculiar velocity of order 1000[km/sec] at this scale. Therefore, at least, the 
quantity T cannot be interpreted as the ordinary temperature as defined from 
the velocity dispersions. One of our next task will be to elucidate the meaning 
of T. 

In relation with the astronomical velocity distributions, the authors[10] 
claim that the velocity distribution of the clusters of galaxies can be well fitted 
by the Tsallis distribution with the parameter q = 0.23^q'q5, which is appar- 
ently different from our negative value for galaxy distributions. Further study 
on the velocity distributions in various scales (galaxies, clusters, super-clusters) 
would reveal the origin and evolution of the large-scale structure of Universe. 
We would like to report these results in the near future. 
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Figure 2: Probability functions / (1) through / (8), from left to right and top to 
bottom. CfA II South Observations (dots with error bars) and our calculations 
(Solid lines) are plotted. 
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